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O ' ABSTRACT 

CO 

' Aims. We discuss and characterize micro-lensing among the 3 brightest lensed images (A-B-C) of the gravitational lens system RXS 
J 1 1 3 1- 123 1 (a quadruply imaged AGN) by means of long slit optical and NIR spectroscopy. Qualitative constraints on the size of different 
emission regions are derived. We also perform a spectroscopic study of two field galaxies located within 1.6 arcmin radius from the lens. 
Methods. We decompose the spectra into their individual emission components using a multi-component fitting approach. A complementary 
decomposition of the spectra enables us to isolate the macro-lensed fraction of the spectra independently of any spectral modelling. 

■ Results. 1. The data support micro-lensing de-amplification of images A & C. Not only is the continuum emission microlensed in those 
Or images but also a fraction of the Broad Line emitting Region (BLR). 

q 2. Micro-lensing of a very broad component of Mg n emission line suggests that the corresponding emission occurs in a region more compact 
}_i than the other components of the emission line. 

3. We find evidence that a large fraction of the Fe n emission arises in the outer parts of the BLR. We also find a very compact emitting region 
OS ■ in the ranges 3080-3540 A and 4630-4800 A that is likely associated with Fe n. 

^ 4. The [Oin] narrow emission line regions are partly spatially resolved. This enables us to put a lower limit of ~ 1 10/7 1 pc on their intrinsic 
. ,— 1 , size. 

■ 5. Analysis of Mg 11 absorption found in the spectra indicates that the absorbing medium is intrinsic to the quasar, has a covering factor of 20%, 
' and is constituted of small clouds homogeneously distributed in front of the continuum and BLRs. 

6. Two neighbour galaxies are detected at redshifts z = 0.10 and z = 0.289. These galaxies are possible members of galaxy groups reported at 
those redshifts. 
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1. Introduction 

RXS Jl 13 155.4-123 155 (hereafter Jl 131) is one of the nearest 
confirmed multiply imaged AGN. The source at z s = 0.66 is 
lensed by an elliptical galaxy at zi = 0.295 (Sluse et al.[2003). 
This system is a long axis quad with an image configuration 
very similar to B1422+231 (Patnaik et al. 1992): three merg- 
ing images (B-A-C; typical of a source lying close to a cusp 
caustic) face the faint saddle-point image D lying close to the 
lensing galaxy G (Fig. [T}. 
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J 1 1 3 1 has several observational advantages compared to 
the other lens systems, namely an integrated magnitude R ~ 
15.3 and a large angular separation between the lensed images 
(A9 ~ 3.6"). These characteristics make J 1 1 3 1 an ideal target 
for spectroscopic studies (e.g. Sheinis 2006) and optical flux 
monitoring (e.g. Morgan et al. 120061 COSMOGRAIL project). 
Additionally, the low redshift of the source quasar enables one 
to observe the Mgn, U/3 and [On] AA4959, 5007 emission lines 
in the optical range. This allows the comparison of the flux 
ratios in three different emitting regions, namely the AGN con- 
tinuum, the broad emission lines (e.g. H/3) and the narrow 
emission lines (e.g. [Om] /L14959, 5007). Because of their dif- 
ferent angular sizes, these three emitting regions will not be 
similarly affected by micro-lensing due to stars in the lensing 
galaxy or milli-lensing by more massive substructures (typi- 
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Fig. 1. Image of Jl 131 (0.8"seeing) with superimposed slit ori- 
entations #1 and #2. Identification of the quasar lensed images 
(A, B, C, D) and of the lensing galaxy (G) are also reported. 



cally M > 10 M ). Therefore, the comparison of the flux ratio 
measured in those three regions may offer a robust diagnostic 
of micro/milli-lensing at work in the system. 

Although milli-lensing had been early suspected in Jl 13 1 
(e.g. Keeton et al. [20031 , Sluse et al. (120061 Paper I) gave 
stronger evidence for the presence of micro-lensing, on the ba- 
sis of a multi-wavelength study of the flux ratios observed at 
different epochs. Even when tentatively corrected for micro- 
lensing, the observed flux ratios were still impossible to repro- 
duce with simple smooth lens models. Claeskens et al. ((2006) 
showed in Paper II that multipole-like models could formally 
solve the anomaly. For this purpose, they constructed lens mod- 
els constrained by the relative astrometry of the QSO images 
provided by the HST images and by substructures present in the 
nearly complete Einstein ring . This also enabled them to recon- 
struct the host spiral galaxy of the source AGN and to study in 
detail its properties. In the present paper, we take advantage of 
long slit spectroscopy of Jl 131 to further study micro-lensing 
in this system and phenomenologically describe its effect on 
different emitting regions and in particular on the broad line 
emitting region (BLR). 

Abajas et al. (2002} have shown, using the Kaspi et al. 
(2000) relation linking the BLR size to the quasar luminosity, 
that the BLR of less luminous quasars is likely small enough 
(i.e. a few micro-arcseconds) to be affected by micro-lensing. 
Because the source in J 1 1 3 1 is rather a (bright) Seyfert than a 
luminous quasar, one may thus expect to observe differential 
microlensing of its BLR. Thanks to the large spectral coverage 
of our data (i.e. 2400 < A < 6900 A ; rest frame) we looked for 
micro-lensing of the BLR and compared its effect on the main 
Balmer lines and on the Mg n line. Additionally, we compared 
the flux ratios in these lines to those in several narrow emis- 
sion lines (NELs), namely [Nev], [On] and [On]. Finally, we 
studied micro-lensing of the blended Fe n emission all along the 
spectra. The location of the Fe n emitting region is still hotly de- 
bated (e.g. Zhang et al. 2006b) and micro-lensing of the latter 



offers the opportunity to compare the size of the Fe n emitting 
region with the BLR size. 

The structure of the paper is the following. We describe the 
observations and the spectrum extraction in Sect. [2] We per- 
form a preliminary analysis to explore the data in Sect. [3] and 
provide more quantitative measurements of the spectral differ- 
ences by means of a multicomponent spectral decomposition 
in Sect. [4] We propose a comprehensive micro-lensing sce- 
nario in Sect. [5] and discuss the implications of the differential 
micro-lensing of the broad emission lines (BELs) in Sect. [6] We 
present some additional results concerning the extended nature 
of the NELs, the absorption systems observed in the QSO spec- 
tra in Sect. [7] together with the redshifts of the lens and of two 
neighbour galaxies. We comment on some working assump- 
tions in Sect.[8]and summarize the results in Sect. [9] 

Unless explicitly stated, we adopt Hq = 70 km s _1 Mpc" 1 , 
fio = 0.3 and Ao = 0.7; all the^- 2 values reported are reduced 
values. 

2. Observational overview 

2.1. Data 

We obtained long slit spectra of Jl 13 1 with the FORS2 in- 
strument mounted on the Kueyen (UT2) ESO Very Large 
Telescope on April 26th 2003. These data consist of two se- 
ries of spectra with the 1" slit oriented along the Jl 131 lensed 
images B-A-C and A-D respectively (hereafter orientations #1 
and #2; see FigQ]). The Longitudinal Atmospheric Dispersion 
Corrector (Avila et al. 1997) has enabled us to keep the slit 
centering within ~ 0.2"all along the wavelength range, re- 
ducing the slit losses in the blue range. The grisms GR600B 
and GR600I+OG590 blocking filter have been used, allow- 
ing to cover the wavelength ranges 3890 < A < 6280 A and 
6760 < A < 88 10 A with a resolving power around 800 and 
1500 at the central wavelength. The CCD camera (SR collima- 
tor) was unbinned (0.063"/pix in the spatial direction) in order 
to ease the spatial deblending of the spectra. Long slit spec- 
tra (with slit orientations identical to FORS spectra) have also 
been obtained in the near infrared range (9810 < A < 11390A) 
with the ISAAC instrument placed at the focus of the VLT UT1 
telescope (Antu), 13 days before the FORS spectra. The low 
resolution grating combined with a 1" slit and a SZ filter leads 
to a resolving power of 550. Table Q] gives a short log of the 
observations. 

2.2. Reduction 

The FORS spectra were bias subtracted and flatfielded using 
IRAF 1 . A set of 5 dome flat fields were used. The wavelength 
dependent structure of the flat fields was corrected by fitting a 
high order cubic spline with the response task. The resulting 
flat field is of sufficient quality to correct the fringing appear- 
ing on the raw data. The wavelength calibration was performed 

1 IRAF is distributed by the National Optical Astronomy 
Observatories (NOAO), which are operated by the AURA, Inc., under 
cooperative agreement with the National Science Foundation (NSF). 
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Table 1. Log of the observations for slit orientations #1 and #2 
(Fig[TJ. Col. 1: Observing date; Col. 2: Instrument; Col. 3: Total 
exposure time; Col. 4: Mean seeing during the observations; 
Col. 5: Mean airmass sec(z). 



Date Instrument Exp seeing sec(z) 

(dd-mm-yy) (s) O 

slit #1 

26-04-03 FORS2 (GR600I) 4 x 240 055 LIT 
26-04-03 FORS2 (GR600B) 4x240 0.4-0.7 1.45 
13-04-03 ISAAC (LR+SZ) 6 x 200 0.45-0.62 1.19 

slit #2 ~ 
26-04-03 FORS2 (GR600I) 3 x 430 053 1.07 
26-04-03 FORS2 (GR600B) 3 x 430 0.6-0.8 1.12 
13-04-03 ISAAC (LR+SZ) 14 x200 0.56-0.89 1.40 




thanks to He-Hg-Cd and He-Ar-Ne calibration arcs. Since im- 
ages B and C are not well centered in the slit, we shifted the ze- 
ropoint of their wavelength calibrated spectra by ~ 0.8 A with 
respect to A. Finally, the sky background was removed by fit- 
ting and subtracting a first order Chebychshev polynomial in 
the spatial direction in two adjacent regions of 180 pixels width 
not illuminated by the object. Since no standard star was ob- 
served during the science night, the sensitivity curve was con- 
structed using the two standard stars GDI 08 (for the GR600B 
grism) and LTT7379 (for the GR600I grism) observed on April 
15th 2003 and May 25th 2003, respectively, through a 2.5" 
width slit. The same reduction procedure was applied to these 
stars, using calibration data associated with these observing 
dates. The stars and the objects were corrected for sky extinc- 
tion with the Paranal extinction estimates available on the ESO 
website. 

For the ISAAC spectra, as is standard practice in the in- 
frared, the object was observed in two positions along the 
slit. The strong and highly variable night sky emission lines 
were effectively removed by subtracting the resulting spectra 
from each other. The 2-D sky-subtracted spectra were then flat- 
fielded, registered, wavelength calibrated and co-added. The 
wavelength calibration was performed using both xenon ar- 
gon arc lines and OH sky lines (Rousselot et al. [2000). The 
slit curvature correction was applied using the fitcoords and 
transform IRAF tasks. Once corrected for the slit curvature, 
the NIR spectra were still inclined with respect to the pixel grid. 
We thus applied an additional rotation of 0.78° to align the 
spectra with the pixel grid and ease the subsequent spectrum 
extraction (Sect. I2.31 l. Finally, the response curve was calcu- 
lated by dividing the reduced spectrum of the B2 type standard 
star HIP88126 by a black body at 20000 K. For both optical 
and NIR data, a noise image was constructed at each step of 
the reduction, using the standard error propagation formula. 

2.3. Spectrum extraction 

Our extraction method is different for slit orientations #1 and 
#2. For frames obtained with slit orientation #1, only the point- 
like images A, B and C were in the slit. We thus fitted three 
Moffat profiles in the spatial direction for each wavelength bin 



Fig. 2. Central panel : Bidimensional Moffat models of the 
spectra of C, A & B in the U/3A4S61 & [O m] AA4959, 5007 re- 
gion. Left and right panels: bidimensional residuals (+3cr) for 
a Moffat and a Gaussian fit. Vertical black lines indicate the 
center of the profiles for C, A, B. 



independently, allowing both the centroid of the Moffat, the 
slope (3 and FWHM to vary. After this preliminary fit of the 
bi-dimensional spectra, we smoothed the fitted centroid posi- 
tions of A-B-C with a 10 A moving average and reprocessed 
the extraction by fixing the positions to the ones found after 
smoothing. The fitted 2D spectra were then compared to the 
observed ones in order to track for extraction artifacts. Figure|2] 
compares the residuals observed in the H [} and [O m] region 
(i.e. 7900 < A < 8350 A) when extraction is performed with 
3 Moffat (left panel) and with 3 Gaussian (right panel) pro- 
files. One can see systematic residuals next to each lensed im- 
age spectrum when Gaussian profiles are used. These residuals 
are indeed observed all along the spectra. On the other hand, 
the residuals are flat at all wavelengths when Moffat profiles 
are used, except around [O m] AAA959, 5007. In order to ex- 
plain the origin of those residuals we constructed a 2D pseudo- 
continuum under [O m] AAA959, 5007 by interpolating the con- 
tinuum level measured just before and after [Om], Once that 
pseudo-continuum had been subtracted from the data, we ob- 
tained 2D spectra of [O m] emission only. The fit of the latter 
with 3 Moffat profiles is significantly worse than the fit of the 
same region for the continuum only (x 1 = 1.01 for the con- 
tinuum and x 1 = 3.71 for [Om] -U4959, 5007). The results are 
even worse if we impose the Moffat parameters for [O m] to be 
equal to those found for the continuum (x 2 - 4.68). This is a 
hint that the observed [O m] AA4959, 5007 emission is not point- 
like but is partially resolved. The complete extracted spectra 
of A, B and C are displayed in Fig. [3] The flux "lost" in the 
[Om] region due to imperfect extraction is only 0.5% of the 
[O in] /U4959, 5007 emission in A+B+C. Therefore, our ex- 
traction method does not introduce significant error on the flux 
measurement in [O m] due to its spatial extension. 

After extraction, we constructed a synthetic optical spec- 
trum (sampled on a pixel grid of 1.33 A) of each lensed im- 
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9800 10000 10200 10400 10600 10800 11000 11200 11400 

Wavelength (Angstroms) 

Fig. 3. Extracted spectra of images A (thick solid red), B (dashed green), C (dashed-dotted blue) and D (thin solid pink) in the 
blue (up), red (middle) and NIR (bottom) ranges. Black dashed line (bottom of the graph) in the blue and red ranges is the 
template of Fe n emission redshifted at z = 0.657 and arbitrarily shifted in flux for legibility. Relative flux ratios between the 
spectra are conserved except for image D. Main emission lines are identified. 



age by joining the blue (GR600B) and red (GR600I) spectra 
on a common grid. Since there is no overlapping domain be- 
tween the blue and red spectra, we assume that the fractional 
loss of flux is similar in each spectrum such that no relative 
flux rescaling between blue and red is necessary. We did not 
attempt to rescale in flux the NIR spectra over the optical ones 
since these spectra were obtained 13 days later. Although A, B 
and C are not perfectly aligned along the slit, we assume that 
the relative flux loss between the lensed images is small. To test 
this, we extracted the flux in the acquisition image along 3 dif- 
ferent rectangular apertures of 0.8, 1.0 and 1.2" widths, each 
oriented like the true slit. We then fitted A, B, C with 3 one- 
dimensional Moffat profiles as we did for the spectrum extrac- 
tion. Following this procedure we found that the measured flux 
ratios are independent of the chosen slit width. Additionally, 
these flux ratios agree within 5% with the ratios of the spectra 
multiplied by the 7?-band transmission curve. 

For the slit orientation #2, we performed classical extrac- 
tions in apertures centered on the lensed images A and D and 
around G, using the IRAF task apall. Unfortunately, the ex- 
tracted spectra of D were significantly contaminated by the 
lensing galaxy while the one of A was likely contaminated by 
images B and C. Although we could adapt the size of the ex- 
traction aperture to reduce the contamination, this precluded a 



correct flux calibration for the spectrum of D. For this reason, 
we do not analyse that image in the remaining. Additionally, 
three other objects are located in the slit. A star (11:31:55.8, 
-12:32:16) and two galaxies gal#\ and #2 (11:31:55.025, - 
12:32:13 and 11:3 1:57.7, -12:32:22) located at 55 and 95" from 
the lens. They were extracted through fixed apertures of 10 
pixel width. These spectra are analysed separately (Sect. I7.3I > 
from those of the lensed images. The latter are described in the 
following two sections. 

3. Exploratory analysis 

The goal of this section is to identify general features in the 
A-B-C lensed image spectra, by means of a phenomenological 
and model independent approach. Although numerical results 
will be obtained, more quantitative results derived from sys- 
tematic fittings will be presented in Sect. [4] 

3.1. Spectrum ratios 

The first insight into the relative differences between the lensed 
spectra is simply provided by computing the ratios of individ- 
ual spectra. Indeed, since gravitational lensing is achromatic, 
one expects the lensed image spectra to be identical up to a 
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On the other hand, the C/B spectral ratio is not flat. First, 
most of the BEL structure has disappeared, with the notice- 
able exception of the Mgn line, in which the flux ratio is 
significantly different from the one in the underlying con- 
tinuum. Second, the flux ratio in the NELs (in particular in 
[O in] AA4959, 5007) is drastically different from the one in the 
continuum (i.e. as measured at the foot of the line "ratio"). 
Finally, an overall, chromatic trend is observed with a flux ra- 
tio increasing from ~ 0.35 in the blue to ~ 0.55 in the red. An 
attempt to provide a coherent interpretation with a minimum 
number of hypotheses is discussed below. 



6000 7000 
Observed wavelength (Angstro 




3.2. The microlensed spectra 
3.2.1. The Fn method 

Assuming the observed spectra F, are simply made of a super- 
position of spectrum Fm which is only macro-lensed and of a 
spectrum both macro and m/cro-lensed, it is easy to ex- 
tract both components Fm and Fm/i by using pairs of observed 
spectra. However, the relative macro amplification M must be 
known. 

Indeed, defining M = Mi/M 2 (> 0) and jj. = /ii/yLi 2 (> 0) 
as the constant macro- and micro-amplification ratios between 
image 1 and image 2, we have: 



6000 7000 
Observed wavelength (Angstroms) 



Fig. 4. Ratio of spectra A & B (top) and C & B (bottom). The 
observed (flat) ratio R = M/i is plotted with a dotted line. R = 
0.63 for the A/B ratio and R = 0.33 for the C/B ratio. 



magnification factor. Therefore the spectral ratio between two 
images is expected to be flat. However, extinction in the lens- 
ing galaxy may differently affect the lensed images, as well 
as micro/milli-lensing produced by substructures in the lens, 
such that chromatic trends may appear in the spectral ratios 
(e.g. Jean & Surdej 119981 Wambsganss & Paczvriski ll991| ). In 
some cases, micro/milli-lensing can also modify the profiles of 
the BELs (e.g. Abajas et al. 2002). The ratio between individ- 
ual spectra thus provides one with a simple diagnostic of the 
presence of differential extinction or micro-lensing. The A/B 
and C/B ratios are presented in Fig. [4] 

One can clearly see the imprint of the QSO emission line 
spectra in the A/B spectral ratio, indicating that the flux ratio is 
different in the continuum (e.g. 7900-7950 A), the BELs (e.g. 
8000-8170 A) and the NELs (e.g. 8240-8320 A). Differences 
in the pseudo-continuum (i.e. the AGN power law contin- 
uum+blended Fen emission) are also apparent (e.g. 5100- 
6200 A). This suggests that "classic" microlensing is probably 
at work by affecting the compact region emitting the contin- 
uum. However, no chromatic dependence is observed, as illus- 
trated by the flat underlying continuum ratio (tentatively indi- 
cated by the dotted line in Fig.|4j at the value R = 0.63) 



Fx 
F 2 



MF M + M/iFm/j 
Fm + Fm/j ■ 



(D 



The latter equations can be rewritten to extract Fm and F^: 



F m - 

Fm/j 



F\ jM-jiF2 

_ F 2 -Fi /M 
1-H 



(2) 



where fi must be chosen to satisfy the positivity constraint 
F M > and F Ml , > 0. 

It is clear from Eq. (O that Fm/j is directly determined from 
the macro-amplification ratio M and the observed spectra Fi 
and Fi, up to a scaling factor. On the other hand, if /u and // are 
two values of the micro-amplification ratio which satisfy the 
positivity constraints and such that // < fi, it is easy to show 
that: 



F M (M') = F M (M) + — -F 



1 



(3) 



It is obvious from the previous equation that adopting the 
smaller value // is equivalent to add a fraction of the mi- 
crolensed spectrum to the macrolensed spectrum obtained with 
[i. Since an emitting region can only be microlensed or not, the 
only possible choice for fi is the maximum value satisfying the 
positivity constraint. 

The previous decomposition is exact when fj. is constant as 
a function of wavelength, i.e. when the micro-lensing ratio is 
achromatic. However microlensing is a function of the source 
size (which is wavelength dependent) and this function can be 
different for different images, so that even the micro-lensing 
ratio fi can be chromatic. Any such wavelength dependence of 
pi would be propagated into the spectra F M and F M ^,- Although 



6 



Sluse D. et al.: Long slit spectroscopy of RXS Jl 13155.4-123155. 



this is not important as long as the method is simply used to 
identify spectral features affected by some micro-lensing (i.e. 
whose flux ratio is not equal to the macro value), we shall show 
that p ^constant is a reasonable choice. 

In the following, we will call the decomposition of the spec- 
tra into Fm-Fm/i, the F p method. 

3.2.2. The A-B pair 

In order to apply the Fp method to the A-B pair of spectra, 
the macro amplification ratios must be fixed. In the remain- 
der of the paper, we use the hypothesis that they are given by 
the observed flux ratios in the [Om] /U4959, 5007 lines. This 
seems a reasonable choice in the sense that the NELs are pro- 
duced by the largest emitting region, thus the least affected by 
microlensing. The fact that the [O m] AA4959, 5007 emission 
is partially resolved, especially between images C and A (see 
Sect. 12. 31 does not bias the [O m] flux estimate. Indeed, we have 
shown that only a very small fraction (i.e. <1%) of the flux in 
[Om] was lost during the extraction (Sect. l2~3l >. 

Under this hypothesis, first we must have that M = 
M A /M B = 2.1 in order to suppress the [O m] AA4959, 
5007 emission lines in the combined spectrum Fa - MFb- 
Second, by associating image A (resp. B) with number 1 (resp. 
2) in Eqs. (0, we find p = PaIpb = 0.3. The extracted spectra 
Fm and Fm^ are shown in Fig. [5] Note that a constant value of 
p is certainly a good approximation, given the flat A/B under- 
lying continuum (Fig. |4j 

The striking result of this decomposition is that, in addition 
to the continuum, the broadest part of the BEL is also affected 
by microlensing. This might indicate that at least part of the 
BLR is very compact. A more thorough discussion is delayed 
to Sect. [5] after a more quantitative analysis is made in Sect. [4] 

On the other hand, as expected, the unaffected spectrum 
Fm is basically flat and consists of the NEL and narrowest 
part of the BEL. Indeed, not only the [O m] AA4959, 5007 lines, 
but also the core of Mg n, H y and H B BELs are prominent. 
Interestingly, blended Fen emission is also clearly observed 
around Mg n and blue/red-ward of H B. 

3.2.3. The B-C pair 

In the case of the B-C pair of spectra, before applying our ex- 
traction method, we must first discuss the chromatic trend ob- 
served in the ratio spectrum (see Fig. 0}. This trend probably 
comes from image C since the A/B ratio is not affected. It has 
four possible origins: i- chromatic microlensing; ii- differential 
extinction along the line-of-sight; iii- host contamination; iv- 
contamination by the lensing galaxy. 

We exclude significant contamination by the lensing galaxy 
that is about 100 times fainter than image C at the lensed 
images position (Paper II). Since none of the other phenom- 
ena taken separately is capable of accounting for the obser- 
vations, we must search for the minimum combination of hy- 
potheses. Indeed, extinction alone is unable to explain the pres- 
ence of Mg ii in the ratio spectrum and is expected to modify 
more strongly the spectrum slope in the blue than in the red. 
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Fig. 6. The smoothed host galaxy spectrum (dotted) responsi- 
ble for the chromatic trend observed in the C/B spectrum ratio, 
compared with a redshifted spectral template of a Sb galaxy 
(solid ). 



Chromatic microlensing is also expected to mainly affect the 
blue part of the spectrum; finally, the host galaxy is expected 
to maximize its contamination in the reddest part of the spec- 
trum of the faintest image, as observed, but it is not supposed 
to produce a differential effect on Mg n. 

Since the host is indeed observed to be quite bright in the 
red and in the NIR (see Papers I & II), we decided to ex- 
plore the combined effect of the host contamination together 
with a wavelength independent microlensing, as modelled in 
Sect. ED] 

The host contamination probably seen in the spectrum ra- 
tio necessarily comes from a constant, resolved contribution 
(the possibly unresolved host contribution undergoes the same 
macro-amplification as the QSO and is spatially too large to be 
microlensed, so that it cannot affect the slope of the spectrum 
ratio). The contamination by such a resolved contribution is ad- 
ditive and, thus, it must first be removed before applying Eqs. 
©. Since, conversely, the dominant part of the resolved flux 
may be associated with the host and the dominant part of the 
unresolved flux may be associated with the QSO continuum, 
we can write with a good approximation: 

F B = Fqso + Fftost 
Fc - RFqso + F Host , 

where R is the intrinsic C/B ratio due to macro- and micro- 
amplification. Fuost is obtained after smoothing out residual 
features in the solution of system[4] 

Now, the exact choice of the R value is somewhat arbitrary. 
If the host contamination in the blue is assumed to be negligi- 
ble, the value R = 0.35 is chosen, as observed in the blue part 
of the spectrum ratio (see Fig. |4j». A break is seen in the result- 
ing host restframe "spectrum" exactly at the expected position 
of the Balmer discontinuity, giving a good hint that we indeed 
deal with the host galaxy (unfortunately the break falls in the 
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3000 4000 5000 6000 7000 8000 9000 10000 

Observed wavelength (Angstroms) 

Fig. 5. Fraction of the spectrum affected (F^, dotted) and unaffected (Fm, solid) by microlensing (arbitrary units) in the A-B 
spectrum pair (see Eq. (01). The A-band telluric absorption is also identified. 



core of Mg n and H f5 BELs). This would confirm that the BLR 
is very compact. However, as we shall discuss in Sect. |5J this 
decomposition has to be compared with the one of the A-B pair 
and interpreted in terms of the micro-lensing of several lensed 
images. 

Before entering into the possible interpretations, we would 
like to further confirm the previous results by deriving the dif- 
ferent flux ratios in a more quantitative way, based on a sys- 
tematic fitting of the different spectral components (i.e. host, 
QSO-continuum, QSO emission lines, ...). 

4. MCD: Multi Component Decomposition analysis 

Examination of our extracted spectra reveals significant emis- 
sion associated with Fen. An important fraction of the flux 
of the host galaxy is also mixed with the flux from the cen- 
tral AGN (Sect. 13.2.31 ). In order to quantitatively disentan- 
gle the different emission components (which might be differ- 
ently affected by microlensing), we used a Multi-Component 
Decomposition (MCD) approach (e.g. Wills, Netzer and 
Wills 1985 1 Dietrich et al. 120031) . This method assumes our rest 
frame spectra to be a superposition of: (1) a power law contin- 
uum (F v ~ v°), (2) a Balmer emission continuum, (3) a pseudo- 
continuum due to the merging of Fe n emission blends, (4) a 
galactic template for the AGN's host galaxy, and (5) an emis- 
sion spectrum due to the other individual broad emission lines. 
The fitting of the spectra is a classical least square minimization 
using a Levenberg-Marquard based algorithm adapted from the 
Numerical recipes routine (Vetterling et al. 1993). We give 



Table 2. Summary of the macro-amplification ratios M, micro- 
amplification jj and continuum ratio R values found in the ex- 
ploratory study. 





A/B 


C/B 


M 


2.1 


1.4 


M 


0.30 


0.25 


R 


0.63 


0.33 



gap between the two wavelength ranges covered by our spec- 
tra). As a refinement, adopting R = 0.33 instead of R — 0.35 
adds a slight host contamination in the blue, so that the final 
contaminating "spectrum" is reminiscent of the spectrum of a 
Sb galaxy (see Fig.|6]l. In the first (resp. second) case, the host 
flux at 5000 A restframe corresponds to ~ 50% (resp. ~ 75%) 
of the QSO flux observed in image C at the same wavelength. 
However, in both cases, the C/B spectrum is flattened and no 
significant difference is observed in the subsequent analysis. 
Note also that correcting a posteriori for the host contamina- 
tion in the A-B pair does not significantly affect the results de- 
rived in Sect. 13.2.21 

Once the B and C spectra are decontaminated (using R = 
0.33), the spectral decomposition is performed in the same way 
as for the A-B pair. Associating image C (resp. B) with number 
1 (resp. 2) in Eqs. ©, we find M = M C /M B = 1.4 and p = 
PcI^b - 0.25; the extracted spectra are shown in Fig. [7] By 
construction, the unaffected spectrum is flat. The BEL appears 
here nearly totally affected by microlensing (except parts of the 
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3000 4000 5000 6000 7000 8000 9000 10000 

Observed wavelength (Angstroms) 

Fig. 7. Fraction of the spectrum affected (Fify, dotted) and unaffected (Fm, solid) by microlensing (arbitrary units) in the B-C 
spectrum pair, after removing the observed host contamination (see Eq. (|2])). The A-band telluric absorption is also identified. 



hereafter the technical details associated with this method. We 
use it to decompose the observed spectra of A-B-C (Sect. |4~21 ) 
and discuss the errors in Sect. 14.31 

4.1. Overview of the method 

We performed the fit iteratively. First we fitted components 
(1), (2) and (4) in three rest-frame windows nearly free of 
Fe ii emission, namely /M2600-2700 A, 2950-3 100 A and 5400- 
5500 A (Natali et al. I1998I I. Secondly, we imposed a null 
weight to the emission line regions and fitted components (1) 
to (4) using the previous results as initial conditions. Finally we 
included the emission lines in the fit. 

4.1 .1 . The quasar continuum 

We modelled the quasar emission as: 

i -(?+<*) 

F A = F x — + Ff c (5) 

where Aq = 5000 A, a is the canonical power law index and 
F^ c is the Balmer continuum intensity. To estimate the Balmer 
continuum from our fits, we have used the empirical distribu- 
tion by Grandi d 1982b which assumes gas clouds of uniform 
temperature T e that are partially optically thick. In this case the 
Balmer continuum spectrum (A < 3646 A) can be described by: 

F* D = F m B A {T e Xl-e^*), A<A BE (6) 



with Bx(T e ) as the Planck function at the electron temperature 
T e , t,\ as the optical depth at A, and F BE as a normalized es- 
timate for the Balmer continuum flux density at the Balmer 
edge A^E = 3646 A. At wavelengths A > Abe higher order 
Balmer lines are merging into a pseudo continuum, yielding 
a small rise to the Balmer edge. We did not attempt to model 
this region and kept F Bac = at 4150 A > A > A BE . In or- 
der to ease the minimization, we fixed the electron tempera- 
ture and the optical depth t,\ to realistic values (i.e. tx - 2.0 
and T e = 15000 K; Wills et al. [19851 Dietrich et al. 120031 . 
The region 3080-3540 A is imperfectly modelled using Eq. ([6}. 
We observe a significant excess of emission w.r.t. the model. 
This one is likely explained by blended Fe n emission for which 
there is presently no available template and which are poorly 
reproduced by photo-ionization models (e.g. Veron-Cetty et al. 
2006). We thus excluded this region from the fitting procedure. 

4.1.2. The Fell emission 

Following many authors (e.g. Wills et al. U985l Laor et al.Q.997 , 
Dietrich et al. 2003, Tsuzuki et al. 2006), we used template 
spectra of Fe n emission instead of theoretical emission models. 
Namely, we fitted the Fe n emission in our spectra by using an 
UV Fen template (Vestergaard & Wilkes 1200 11 1 and an optical 
Fe ii template (Veron-Cetty et al. 2004). We also tested the opti- 
cal template of Boroson and Green (119921 1 but this did not mod- 
ify significantly our results. We prefered the optical template 
of Veron-Cetty et al. especially due to its larger spectral cover- 
age. Since UV photons are generally absorbed and re-emitted 
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at optical wavelengths (for high optical depth), optical and UV 
Fen emission are often anti-correlated (Joly M., private com- 
munication). For this reason, we have not fixed the intensity ra- 
tio between the UV and optical Fen templates. The templates 
have been broadened in the log space (in order to keep the ve- 
locity width constant over the whole wavelength domain) using 
the FWHM of the H /3 emission line. The broadened templates 
are shown in Fig. [3] 

For technical reasons, in the Fe nuv template of Vestergaard 
& Wilkes, the Fen level beneath the central part of the 
Mgn emission (i.e. 2787 < A < 2802 A) is 0. We have thus 
constructed a modified Fen template for which we consider 
a non zero Fen level between 2752 < A < 2832 A. We esti- 
mate this level by interpolating in the template the Fe h level 
between the two closest maxima enclosing the zero Fen hole. 
These maxima are found at 2752 and 2832 A. Both templates 
have been used during our decomposition. The effect of this 
modified template on the results is discussed in Sect. 14.2.21 

4.1 .3. Host galaxy 

We used a Sb galactic template (Kinney 1993) to model the 
host galaxy. The latter choice is motivated by the detailed study 
of the host performed in Paper II. It is also supported by the 
empirical host spectrum we have retrieved in Sect.[3]and which 
has a continuum emission compatible with a Sb type galaxy 
(Fig.il. 

4.1 .4. Emission lines 

During the last iteration of our fitting procedure, the emission 
lines were fitted simultaneously with the other components as 
a sum of Gaussian profiles. A minimum of two profiles was 
considered for the main emission lines. More details are given 
in Sect. l432l 

4.2. Application to A-B-C spectra 

We present in this section the results of the MCD method ap- 
plied to the spectra of images A, B and C in order to estimate 
the flux ratios in the different emission regions. Unfortunately, 
we found that there is a degeneracy between three quantities: 
the host galaxy relative flux, the power law index a and the 
Balmer continuum level F^ ac . Due to the large range of ac- 
ceptable values for a implied by the degeneracy, we could not 
disentangle the host galaxy emission from the QSO emission 
and we could not assess the presence (absence) of chromatic 
variation of the quasar continuum (due either to microlensing 
or differential extinction). In the following, we thus decide to 
make a prior on a in order to break the degeneracy. This has no 
major effects on our conclusions as explained in Sect. I4.3I 

Motivated by the results of Sect. [3] we imposed a to be 
identical for each image. We choose a = -0. 125 which is at the 
middle of the overlapping range of acceptable values for a in A, 
B, C, when no prior is imposed. This value is also similar to the 
best value found for image B. This value is in good agreement 
with the mean power law index a = -0.33 + 0.59 found by 



Natali et al. (1 19981 1 in the range 1200-5500 A or with the power 
law index of the composite SDSS quasar spectra measured in 
the range 1300-5000 A (i.e. a = -0.44; Vanden Berk et al. 
120011 . Results are given in Table[3] 

4.2.1 . The continuum 

The flux ratios between the power law continua are directly 
obtained from the MCD results (Table [3]) and are reported in 
Tabled We find F A /F B = 0.55+0.02 and F C /F B = 0.31+0.01, 
in good agreement with the rough measurements performed on 
the direct spectral ratios (Fa/Fb = 0.63 and Fc/Fb = 0.33, 
see Tabled Since the estimate of the host galaxy contribution 
to the observed spectra is very different with both methods (i.e. 
180% of the QSO flux in C from the MCD method and 75% 
from the direct ratio), we are confident that the error on the 
host galaxy level does not dramatically bias our continua flux 
ratios 2 . 

4.2.2. The emission lines 

Faint NELs were fitted with a single Gaussian 3 profile while the 
bright [Om] AA4959, 5007 were fitted with 2 components. The 
latter decomposition enables us to take into account the line 
asymmetry of the [O m] AA4959, 5007 emission lines. Since 
the FWHM of the narrowest component is close to the spectral 
resolution, one should not over-interpret this decomposition. In 
the following, we will not use the flux measured in the individ- 
ual components but rather the total flux in the line. 

Regarding the BELs, as it is a common practice, we decom- 
posed each line into two broad Gaussian profiles (e.g. Veron- 
Cetty & Veron|2000] Sulentic et al. l2000l Romano et al. [l996l > 
and we added a narrow component when necessary. We en- 
forced the FWHM and A c of the various sub-profiles to be 
identical for each lensed image. Motivated by the claimed sim- 
ilarity between the Mg n and Balmer line profiles (e.g. Mc Lure 
& Jarvis 2002 ), we enforced the velocity width of the 2 broad 
components (i.e. FWHM ~ 2000 and 4300km s" 1 ) to be iden- 
tical -within 50 km s _1 - for each broad line. However, like 
e.g. Corbett et al (2003), we found necessary to add a very 
broad component (FWHM- 22 000km s~') to the Mgn pro- 
file. We did not find any evidence for such a component in the 
Balmer line profile. We supposed that the different decompo- 
sition found for Mg u could be due to the hole of Fe n beneath 
Mg ii in the Fe uuv template. However, the use of a modified 
Fe u template in which we "filled" the Fe n hole (Sect. I4.1.2I I 
does not prevent the presence of a very broad component in 
the Mg n profile. The fit of Mg n was even found to be worse 
with this modified template and especially, the peak of the 
Mg ii emission could not be fitted correctly anymore. Finally, 
we modelled the faint Mg n narrow absorption line doublet seen 

2 Note that it neither affects the Fy decomposition of the spectra. 
Namely, a similar Fyi decomposition is obtained with these two differ- 
ent host estimates. 

3 Due to the low contrast between [Neiv] /12419 and 
Hy+[Oin]/M4340, 4363 emissions w.r.t. their underlying continuum, 
the fit of those lines is inaccurate and not presented. 
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Table 3. Multi-component best fit of the spectra when a is forced to be identical in A, B, C. Col. 1 is the image ID, Col. 2 is the 
power law index a of the quasar continuum, Col. 3 gives the flux level of the power law continuum at A — 5000 A (rest-frame), 
Col. 4 gives the flux level of the Balmer continuum, Col. 5 gives the relative flux of the host galaxy at 5000 A normalized by 
F™ h (in B), Col. 6 is the level of the optical template of Fe n, Col. 7 is the level of the UV template of Fe n. All the quantities 
(except a) are normalized to image B. 



image 


a 


c-PWL 
r 


/?Bac 


-Fhost.O 


Fe 


"opt 


Feii uv 


A 


-0.125 


0.55± 0.02 


0.94+ 0.09 


0.58± 


0.02 


1.91=1 


: 0.03 


2.29± 0.08 


B 


-0.125 


1.00± 0.01 


1.00± 0.09 


0.77± 


0.02 


1.00:1 


z 0.03 


1.00± 0.08 


C 


-0.125 


0.31± 0.01 


0.41± 0.09 


0.54± 


0.01 


0.43d 


: 0.03 


0.82+ 0.06 



Table 4. Multi-component fit of the main emission lines normalized to the total flux in Wj3 for image A. Note that Ha normal- 
ization is arbitrary due to the wrong flux scaling of the NIR spectra of Jl 131 w.r.t. the optical spectra. 



ID 


A c 


FWHM 


F A 




Fc 


z 


Mg iiA2798 


BC1 


4636.4 


2050 


0.281+0.006 


0.077±0.006 


0.081+0.004 


0.657 




BC2 


4638.4 


4320 


0.867±0.009 


0.557±0.008 


0.251+0.006 


0.658 




VBC 


4629.4 


21830 


0.728+0.013 


0.646±0.009 


0.215±0.006 


0.655 




TOT 






1.842+0.017 


1.253±0.014 


0.534±0.009 






no VBC 






1.148±0.010 


0.634±0.010 


0.332±0.007 




H£ ,14861 


NCI 


8039.6 


460 


0.040±0.001 


0.012±0.001 


0.016±0.001 


0.654 




BC1 


8054.2 


2000 


0.345±0.005 


0.134±0.005 


0.091±0.004 


0.657 




BC2 


8055.4 


4310 


0.615±0.006 


0.477±0.006 


0.159±0.005 


0.657 




TOT 






1.000±0.008 


0.623±0.008 


0.266±0.006 






BC 






0.960±0.008 


0.611±0.008 


0.250±0.006 




Ho- ,16562 


NCI 


10850.8 


550 


0.496±0.015 


0.242±0.015 


0.276±0.015 


0.653 




NC2 


10831.7 


550 


0.591±0.013 


0.244+0.013 


0.293+0.013 


0.654 




NC3 


10884.8 


550 


0.199±0.018 


0.158±0.018 


0.337+0.018 


0.653 




BC1 


10874.5 


2070 


3.595±0.063 


1.686±0.063 


1.122±0.063 


0.657 




BC2 


10872.5 


4300 


5.480±0.078 


4.227±0.079 


1.533+0.080 


0.657 




TOT 






10.361±0.104 


6.557±0.104 


3.562±0.105 






BC 






9.075±0.101 


5.913±0.101 


2.655±0.102 




[O m] ,1,14959, 5007 


NCI 


8202.1 


240 


0.052±0.001 


0.031±0.001 


0.038±0.001 


0.654 




NC2 


8198.8 


660 


0.073±0.002 


0.034±0.002 


0.049±0.002 


0.653 




NC3 


8280.7 


310 


0.214±0.002 


0.110±0.001 


0.150±0.002 


0.654 




NC4 


8275.6 


800 


0.172±0.003 


0.084±0.002 


0.108±0.002 


0.653 




TOT 






0.5 11 ±0.004 


0.259±0.003 


0.345±0.003 




[Ne v] ,13425 


NCI 


5664.3 


520 


0.033±0.002 


0.017±0.002 


0.019±0.001 


0.653 


[O a] ,13727 


NCI 


6171.6 


620 


0.023±0.001 


0.013±0.001 


0.023+0.001 


0.656 



in the Mg n broad emission with 2 Gaussian profiles separated 
by 7.2 A rest frame (Sect.lT2l. 

For images B and C, our best fitted profile of the broad 
emission lines was indistinguishable (in terms of Ax 2 ) from 
the profile obtained when no prior was imposed on the decom- 
position. For image A, when no prior was imposed on the fit, 
we found that the best decomposition of the Balmer lines was 
made of two broad profiles separated by typically 1000 km s _1 
suggesting some asymmetry in the line profile. 

We report in Table [4] the results of the fits performed with 
the MCD method. Figure [8] shows the continuum subtracted 
emission lines in each lensed image, together with their mod- 
elled profile. 

The flux ratios calculated using the results of Table [4] are 
similar for Ha,H/3 and Mg □ BELs. The flux ratios calculated 
for the narrow [On] AA4959, 5007 and [Nev] ,13425 emission 
lines agree within the error bars but are different from those 
obtained in the BELs. This independently confirms our previ- 



ous finding that both Fa/Fb and Fc/Fb have different values 
for the BELs and for the NELs. Note that due to the excel- 
lent agreement between the flux ratios measured in [O m]/l4959 
(NC1+NC2) and [Om]^5007 (NC3+NC4), we do not distin- 
guish between these two lines in the following and rather con- 
sider the total flux [Om] tot measured in the [Om] doublet. 

The most precise flux ratio estimates are obtained for the 
broad component BC of H/3 and for [Oin] t ot- They represent 
characteristic values of the BEL and the NEL respectively and 
are reported in Table [5] We note that the [O m] flux ratio values 
(F A /F B ([Om]) = 1.97 and F c /F B ([0 m]) = 1.33) found here 
are in good agreement with the macro-amplification ratios de- 
rived in Sect. [3] to match the [Om] emission in each spectrum 
(see Table|2]i. 

Interestingly, for each BEL, the flux ratios found for the 
broadest components (e.g. BC2 or VBC) are different from 
those of the narrowest broad component (BC1) except for 
FaIF c . 
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8150 8200 8250 8300 8350 10650 10700 10750 10800 10850 10900 10950 11000 11050 



Observed wavelength (Angstroms) Observed wavelength (Angstroms) 

Fig. 8. Comparison of the main NEL (left panel) and BEL (right) for A (solid red), B (solid green) and C (solid blue) with their 
multi-component fitted model (superimposed with dotted lines). Hfj, Mgu and [O m] /L14959, 5007 are continuum subtracted 
profiles. For H a, [Ne v] and [O n] , we have normalized the apparent pseudo continuum to the one of image C. 



The Fe u opt flux ratio (Table [3]l is similar to the Hfj ratio for 
the image pair B-C and to the [O m] ratio for the image pair A- 
B, as can be easily seen in Table [5] There is an indication for 
a similar behaviour in Fe mjv but the flux ratios deviate more 
significantly from those obtained in Hfj and [Om]. 

Finally, a systemic redshift z S ystemic = 0.657 + 0.001 has 
been determined from the Mg h emission. The measured nar- 
row emission lines are all blue-shifted w.r.t. the systemic red- 
shift (znel = 0.654 + 0.001) except [On] emission (Z[oii\ - 
0.656 + 0.001). We also notice that the red wing of the emission 
line profile of [O n] is significantly more pronounced in image 
A (Fig. [8]), probably because of some underlying Fen emis- 
sion (which may deviate from the template we used) and/or to 
emission associated with star formation located in a direction 
orthogonal to the slit. 

In summary, the MCD analysis performed in this section 
confirms the main results unveiled in Sect. [3] First, we mea- 
sure significantly different flux ratios for the continuum emis- 
sion, the BELs and the NELs. Secondly, we confirm differ- 
ential magnification of the broadest component of the BELs 
compared to the core of the emission lines as shown in e.g. 
Fig. [5] Finally, we find that the flux ratio Fc/Fg (resp. Fa/Fb) 
in Fe n is similar to the one of the BEL (resp. NEL). 



Table 5. Summary of the flux ratio values measured with the 
MCD method for several important spectral features. 





A/B 


C/B 


[Om] 


1.97±0.03 


1.33+0.02 




1.57±0.06 


0.41+0.01 


Fen opt 


1.91+0.08 


0.43+0.04 


Continuum 


0.55+0.02 


0.31+0.01 



4.3. Error analysis 

The MCD method provides us with lcr formal errors calculated 
from the covariance matrix in the Levenberg-Marquard min- 
imization routine (Vettering et al. I19931 I. The main limitation 
of the MCD technique comes from a degeneracy between the 
power law index a (and the level F™ h ), the Balmer continuum 
emission (F Bac ) and the galactic template contribution. This 
means that these quantities cannot be retrieved unambiguously. 
Before fixing the same value for a in each image, we have esti- 
mated the range of acceptable values for a and F™ h by fixing 
the host galaxy level to two extreme values beyond which the 
fit becomes unacceptable. For that range of acceptable contin- 
uum spectral decompositions, we found that the fit of the Mg n, 
H fj, [O m] /1/14959, 5007 and Fe n op t remained unchanged. On 
the other hand, the Fenuv intensity was retrieved with only 
20% accuracy. Consequently, imposing the same value for a 
in each lensed image is safe and does not induce systematic 
error on the fit of the emission lines (including Fe n). 

However, we want to be more cautious regarding the re- 
sults of the [Nev]/13425, [On]/l3727 emission lines. Indeed, 
we have no satisfying model for the pseudo continuum un- 
der these lines. Therefore, the Gaussian fitting of these lines is 
performed independently using a flat pseudo-continuum under 
the line. Because of the likely presence of Fe n in that pseudo- 
continuum, our procedure may introduce a systematic bias on 
the fitted line intensity. 



5. Towards a microlensing scenario 

In this section we look whether the difference in flux ratios ob- 
served in the NELs, the BELs and continuum (Sect. [3] & I4.21 > 
can be understood by means of micro-lensing occurring in one 
or several lensed images. We first identify a possible scenario 
following a quantitative approach and then discuss the implica- 
tions. 
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Table 6. Micro-amplification ratio Hi/fij between images i 
and j calculated by normalizing the flux ratio measured 
for a given emission region by the flux ratio measured in 
[O in] AA4959, 5007 (Sect. [5J. For completeness, we have also 
added the micro-amplification ratios computed from optical 
and X-ray imaging at different epochs (mm/yy). Note that 
FWHM(BCl) ~ 2050km s _1 ; FWHM(BC2) ~ 4300km s -1 ; 
FWHM(VBC) ~ 22 000km s -1 . 



ID 


-. 

i^B 


j 


j 

f^C/^B 


[W nij /1/147J7, juu / t t 


1 nn+ n n? 

1 .UUIt U.Ui 


i nn+ n rn 

l.UVH U.UJ 




[Ne v]/l3425 


0.97± 0.10 


1.16+ 0.13 


0.83± 0.11 


[On]-l3727 


0.88± 0.11 


0.68± 0.08 


1.29+ 0.12 


Mgn/12798 BC1 


1.85+ 0.14 


2.33± 0.21 


0.79± 0.09 


BC2 


0.79± 0.02 


2.34± 0.08 


0.34± 0.03 


VBC 


0.57± 0.02 


2.29+ 0.09 


0.25± 0.04 


TOT 


0.75± 0.02 


2.33± 0.07 


0.32+ 0.03 


no VBC 


0.92+ 0.02 


2.34± 0.08 


0.39± 0.03 


H/U4861 NCI 


1.72± 0.21 


1.66± 0.23 


1.03± 0.14 


BC1 


1.31+0.05 


2.56± 0.15 


0.51+0.06 


BC2 


0.65± 0.01 


2.61± 0.10 


0.25± 0.04 


TOT 


0.82± 0.02 


2.54+ 0.08 


0.32+ 0.03 


BC 


0.80± 0.02 


2.59± 0.09 


0.31± 0.03 


Fen pt 


0.97± 0.04 


3.04± 0.28 


0.32± 0.10 


Fen uv 


1.16± 0.12 


1.89+ 0.18 


0.61± 0.12 


Continuum 


0.28± 0.01 


1.23± 0.07 


0.23± 0.04 


R band 04/03 


0.31 


1.24 


0.26 


R band 04/04 a 


0.46 


1.69 


0.27 


R band 04/05 b 


0.58 


1.43 


0.41 


X Ray 04/04 c 


0.09 


0.45 


0.20 



" Paper I 

* Morgan et al. ( 120061 
' Pooley et al. ( 120061 



5.1. A possible microlensing scenario 

Because the narrow line emission takes place in a region larger 
than the BLR and the continuum, it is likely to be much less 
affected by micro-lensing. Therefore, as already mentioned in 
Sect. [3] we can reasonably assume that the macro-amplification 
ratios are close to the flux ratios 4 in [Om] AA4959, 5007. So, we 
can write: 



M: 



F, 



([Om]), 



(7) 



where M, {Mj) is the macro-amplification of image i {j) and 
Fj/Fj is the flux ratio in [O m] calculated for images i and j. On 
the other hand, if micro-lensing /v, (jij) occurs in image i(J), we 
can write for a given emission region: 



F ir ... m Mi 

— (emission) = . 

Fj VjMj 



(8) 



Consequently, if we normalize the flux ratio estimated for a 
given emission region by the flux ratio in [Om] AA4-959, 5007, 



4 The flux ratios measured in [Om]/U4959, 5007(F A /F B = 1.97, 
Fc/Fb = 1.33) slightly differ from those obtained with the SIE+y 
fiducial model used in Paper I (F A /F B = 1.65, F C /F B = 0.9) indicat- 
ing some inadequation between the model and the data (see Sect. l8.lt . 



we infer the micro-amplification ratio fii/fii for that region. 
Table [6] displays the micro-amplification ratio //,■//// for the 
main emission regions computed using the results of Table [4] 
(we have not reported the results for H a because they are sim- 
ilar to those obtained for HB). A first quick look at Table [6] 
shows that, except for the NELs (namely [On], [Nev]), the 
Hillij ratio is not compatible with 1 for any value of (/, j). This 
confirms what we had already revealed in Sect. [3] namely that 
micro-lensing is occurring in more than one image. 

We now examine whether the results are compatible with 
microlensing occurring in two images. A priori, there are three 
possible scenarios, depending on whether the un-affected im- 
age is assumed to be A, B or C. Thus, each scenario corre- 
sponds lo ha - 1, fJ-B = 1 ot fie = 1, respectively. From Table|6] 
it is then easy to derive the respective micro-amplification val- 
ues (/j.b,iuc), (fiA,Hc) an d (jj.a,Hb), both in the continuum and 
in the broad components BC of HB (see values in bold face 
in Table [6}. A robust rule to reject a scenario is the following: 
the BLR being spatially larger than the continuum emitting re- 
gion, it cannot be more affected by microlensing than the latter. 
In other words, the micro-amplification value must be closer 
to 1 in the BLR for the two affected images in the considered 
scenario. It is straightforward to check that the above rule is 
satisfied only in the case when image B is not affected. We find 
in that case that (ha,iic) - (0.28, 0.23) in the continuum and 
(fiA,Hc) = (0.80, 0.31) in HB, which means a de-amplification 
of both images A and C. Both the continuum emitting region 
and the BLR are affected. We now discuss these results. 

5.2. Micro-lensing of the BLR 

As already stated in the introduction (Sect. Q]), microlensing 
of the BLR is not surprising in the case of J 1 131. Here we 
demonstrate that microlensing has indeed the correct scale to 
affect the BLR. The Einstein radius Re of a star of mass M is 
given by: 



Rr 



4-GM D h D os 



D„ 



= 14.3 



Mh~ 



It - days, 



(9) 



where D os , D[ s , D a \ are the angular-size distances between ob- 
server and source (os), lens and source (Is) and observer and 
lens (ol). 

On the other hand, Kaspi et al. (2005} have deduced a 
relation between the size of the BLR and the luminosity 
L,i(5100 A) of the QSO. An improved relation correcting from 
the host galaxy contamination has been derived in Bentz et al. 
(2006). Using the magnitude of image B in the F814W filter 
measured with HST (Paper II), and assuming an amplification 
factor of 10 for that image, we can estimate the likely size of the 
BLR using Bentz et al.'s relation. We find that the luminosity 
AL A (49QQA) = 7.91 10 43 /T 2 erg/s which translates into a likely 
size for the H B emission line region of ~ 40-50 lt-days (using 
h = 0.7; 20-40 lt-days using Kaspi et al.'s relation). Although 
this can only be a rough estimate due to the bias on the lumi- 
nosity introduced by the presence of [O m] + HB emission in 
the F814W filter and due to the uncertainty on the macro-lens 
amplification, we see that a few solar mass micro-lens has an 
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Einstein radius large enough to affect a significant fraction of 
the broad emission line region. 



5.3. Micro-de-amplification of image A 

The micro-de-amplification we find for image A is an indepen- 
dent confirmation of the results found on the basis of flux ra- 
tio temporal variability in the optical (Paper I). They also fully 
support X-ray observations. Indeed, because there is evidence 
that the QSO X-ray emitting region is smaller than the contin- 
uum emitting region (e.g. Poole v 120061 1. the de-amplification 
is expected to be stronger in the X-ray than in the optical, at 
a given epoch. This was clearly the case in April 2004, with 
Pa,x = 0.09 and pAjt = 0.46, as reported in Table [6] from 
data in Pooley et al. (2006) and in Paper I, respectively. Finally, 
the parity of image A being negative (saddle point image), de- 
amplification can easily be produced by substructure(s) of any 
mass (e.g. Keeton 2003). Given the time variability detected 
in Paper I, the microlensing regime is favoured rather than the 
milli-lensing one. 



5.4. Micro-de-amplification of image C 

Like for image A, the de-amplification of image C is possibly 
supported by the even stronger de-amplification observed in the 
X-ray domain in April 2004 (i.e. pc.x = 0.20 and pc,R = 0.27, 
see Tabled. Note however that while the X-ray data alone can- 
not discriminate between the micro-de-amplification of C and 
the micro-amplification of B, our data only support the former 
scenario. A sharp drop in the flux of image C has also been ob- 
served between 2004 and 2006 in the X-ray domain (Kochanek 
et al. 2006). Of course, this is not evidence for de-amplification 
in 2003 but simply another sign that C can undergo variable 
microlensing. 

Unlike the parity of image A, the parity of image C is pos- 
itive (minimum image). Consequently, it can only be amplified 
by an isolated substructure. This is not only the case with an 
isothermal clump (Keeton 2003) but also with a point-mass 
perturber in the scenario of micro-lensing (i.e. in the Chang- 
Refsdal 1984 lens model). This is so because for a minimum 
image and k < 1, we always have k + y < 1, where k and y are 
the macro model dimensionless surface mass density and shear 
values at the position of the considered image. Adding a sec- 
ond point mass perturber does not change the result (Grieger et 
al. 1198 9). However, if a non-negligible fraction of the surface 
mass density of the lens is made of micro-deflectors, the prob- 
ability to de-amplify a minimum macro-image increases a lot, 
as illustrated in Fig. 3 of Schechter and Wambsganss (2002). 
Their M10 model (macro-amplification of 10 of a minimum 
image) represents a situation quite similar to the one of the C 
image of J 1 1 3 1 . The observed de-amplification of image C by a 
factor ~ 4 would constrain the fractional contribution of micro- 
lenses to be larger than 30% of the total surface mass density 
of the lens around image C. 



5.5. Conclusion 

We have shown that the most simple micro-lensing scenario 
consistent with the published data on J 1 1 3 1 implies that images 
A and C are both de-amplified by a micro-lens. The latter also 
affects the BELs. We have presented supporting arguments of 
this scenario. 

6. Microlensing as a probe of the QSO structure 

In the previous section, we presented a scenario accounting 
for the discordant flux ratios measured in the continuum, the 
BELs and the NELs. We have shown that the data are com- 
patible with microlensing de-amplification of images A and C. 
Interestingly, both our exploratory and multi-component anal- 
yses (Sects. [3l& 14.2b show that the BEL is nearly completely 
micro-lensed in image C and partially micro-lensed in image A. 
This difference of behaviour is likely due to the larger Einstein 
radius Re of the micro-lens affecting C. Indeed, the micro- 
lensed fraction of a given emitting region increases with the Re 
of the micro-lens (Eq. [9]). In this context, an emission region 
that is micro-lensed in image A should also be micro-lensed in 
image C and should be very compact, while a larger emitting 
region might be micro-lensed only in image C. Finally, even 
larger regions should not be micro-lensed either in A or in C. 
Consequently, we are able to infer information about the BEL 
structure by simply identifying and characterizing the emitting 
regions microlensed in A & C. 

This scenario further supports our assumption that there 
is no chromatic micro-lensing of the optical continuum. 
Indeed, chromatic micro-lensing is expected to be stronger for 
smaller Einstein radii as well as during caustic crossing events 
(Wambsganss & Paczyhski 1991). Due to the smaller Re in 
image A than in C and to the absence of time flux variations of 
the B/C ratio between 2002 and 2004 (Paper I) -indicating that 
there was no caustic crossing event for image C-, we would 
expect chromatic micro-lensing of the optical continuum to 
be stronger in image A than in C. However, the absence of a 
chromatic trend in the A/B spectral ratio (Fig. |4j suggests that 
this effect is negligible. It is thus meaningful to phenomeno- 
logically probe the BLR structure using the MCD analysis of 
A-B-C (Sect. 14. 2t and the results of the Fp decomposition 
(Sect. I3.21 i. The adopted micro-lensing scenario provides us 
with a simpler interpretation of the Fp decomposition of the 
A-B and B-C pairs of spectra. Indeed, the absence of micro- 
lensing in image B means that the Fp decomposition of resp. 
the A-B and B-C pairs unveils the micro-lensed regions in resp. 
A and C. 

We discuss in the next two sections the structural infor- 
mation inferred from the microlensing of the emission line re- 
gions (Sect. 16.1b and of the pseudo-continuum emitting region 
(Sect. ED. 

6.1. Microlensing of the emission line regions 

A first interesting result relates to the differences existing be- 
tween the Mg ii and Balmer emission lines. Indeed, the assumed 
decomposition of the Balmer lines with 2 broad components 
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(BC1 with FWHM ~ 2000 km s _1 and BC2 with FWHM ~ 
4300km s -1 ) did not enable us to reproduce the Mgn profile. 
The foot of that line is clearly larger than the one of the Balmer 
lines and is well modelled with a very broad component VBC 
(FWHM ~ 22 000km s" 1 ). Additionally, we found with the 
MCD method that VBC is significantly de-amplified in both A 
and C (/^a(VBC) = 0.57; yUc(VBC) = 0.25; Table ®. This is 
confirmed by the finding of a very broad component in Fm^(A) 
and Fm)i(C) at the wavelength of Mg n (Figs.[5]and[7j. Because 
it is microlensed in both A & C, the very broad component 
VBC in Mg ii should be emitted in a very compact region. 

We also found systematic evidence that the 2 broad com- 
ponents BC1 and BC2 used to decompose the BELs are dif- 
ferently micro-lensed in image C. We measured yUc(BC2) < 
yUc(BCl) < 1 for Mgn, H/3 and Ha. This demonstrates that 
the broadest component BC2 is more strongly de-amplified 
than the narrowest broad-line BC1, indicating that BC2 is more 
compact than BC1. Additionally, we found with MCD that 
yUc(NCl) = 1 for the narrow H/3 emission. Consistently, this 
narrow component at z — 0.654 is seen in Fm{C) for both 
H/3 and Hy. However, the H/3 and Hy profiles in Fm{C) can- 
not be fitted by a single narrow Gaussian profile at z = 0.654. 
Especially, the foot of the line is too large, likely due to the 
differential micro-lensing of the broad component of the line. 
Contrary to what is observed for the Balmer lines, there is no 
narrow (not micro-lensed) Mgn emission at z = 0.654, but 
a narrow absorption doublet (Sect. 17.21 ). On the other hand, 
a symmetric Mgn emission line is observed at the systemic 
redshift in Fm(C). A Gaussian fit of this profile shows that 
its FWHM is ~ 2500km s . However, we cannot assess 
whether this component is kinematically isolated (i.e. "nar- 
row" Mgn emission doublet with FWHM~ 1000km s _1 ) or 
not (i.e. associated with differential microlensing of the broad 
Mg n line). 

A more speculative result concerns the microlensing of the 
Balmer lines in image A. We observe with the Fp decomposi- 
tion that only the wings of the H/3 and Hy profile are affected 
by some micro-lensing (Fig. 0. On the other hand, the MCD 
analysis reveals that //^(BCl) > 1 while pa(BC2) < 1. This 
apparent difference of micro-lensing regime between BC1 and 
BC2 (i.e. amplification of BC1 and de-amplification of BC2) 
is likely not a real effect but is due to the MCD method which 
assumes that Gaussian profiles describing BC1 and BC2 are 
microlensed as a whole 5 . Such a measurement is indeed com- 
patible with the central part of the line profile being less micro- 
lensed than the wings, such that an excess of flux is measured in 
the core of the line (i.e. BC1), mimicking a micro-amplification 
of the latter. Although only the wings of H y are observed in 
FMfi(A), the profile is not an exact replica of H/3. Especially, 
the reddest fraction of H y is narrower than its homologue in 
H/3. Additionally, it coincides with possible [Om] ,14363 emis- 



3 This assumption is not valid for image A as suggested a priori 
by the evidence for a partial micro-lensing of the BELs in that image 
while micro-lensing de-amplify nearly completely the BELs in C. 



sion at z = 0.654. This would reveal that at least a fraction 6 
of [O m] ,14363 is coming from a region more compact than 
the other narrow emission lines (e.g. H/3, [Om] AA4959, 5007, 
[Nev]) and even as compact as the inner part of the BLR. 
All these results underline that differential micro-lensing of 
the BEL is important in image A but cannot be fully under- 
stood with our phenomenological approach. Nevertheless, we 
interestingly notice that micro-lensing of only the wings of the 
line profile is predicted in Abajas et al. d20021 > for an outflow- 
ing broad emission region. On the other hand, the existence 
of a compact [O m] ,14363 emission is supported by the work 
of Nagao et al. (2001a) who show that a significant fraction 
of [O m] /14363 is likely formed closer to the continuum than 
the other main NELs. Because the critical density of [Nev] is 
only two times smaller than for [O m] ,14363, one would expect 
[Ne v] to be formed in a region similar to [O m] ,14363 (Nagao 
et al. 2001b) and thus be partially micro-lensed, which is not 
observed. 

In summary, we have brought several pieces of evidence 
that the Balmer and Mgn BLRs are different. Especially, we 
have found that there is a very broad compact component in 
Mg ii, not present in the Balmer lines. The narrow emission 
found in H/3 and H y at the same redshift as the other NELs 
(i.e. [Om] ,1,14959, 5007, [Nev]) is absent for Mgn where in- 
stead narrow absorption lines are observed. The micro-lensing 
of the BELs also confirms that the size of the emission region 
is anti-correlated with the FWHM of the line. Finally, it seems 
that the micro-lensing of the BELs in image A resolves the 
BLR in velocity. However, in order to discriminate between 
different BLR emission models, multi-epoch data as well as a 
realistic modelling of the BLR (e.g. Abajas et al. 120021 Lewis 
& Ibata 2004) are necessary. This will be investigated in a fu- 
ture paper. 

6.2. Micro-lensing of the pseudo-continuum 

The Fp. decomposition of the spectra (Fig. [5]&|7]i shows that 
a significant fraction of the Fe n is microlensed in image C but 
not in image A. This is confirmed by the MCD analysis for 
which we measure similar flux ratios in the Fe n opt (and possi- 
bly Fenuv) and in the BELs. This indicates that Fen is emit- 
ted in a region similar to the BLR. This is not surprising since 
Fe n emission lines are generally broadened like the broadest 
part of Balmer lines (e.g. Sulentic et al. [2000). Although our 
result does not prohibit that Fe n emission arises in an accre- 
tion disk as generally believed (e.g. Collin-Souffrin et al. 1980 
Zhang et al. 2006a), it does not favour a scenario where all the 
Fe n is emitted in the very inner part of the disk. Nevertheless, 
we will show hereafter that we have a hint that a fraction of the 
Fe n emission is emitted in a more compact region. 

We can see in Fm^A) that there is an excess of flux just 
blueward of the H/3 emission (i.e. rest frame range 4630- 
4800 A). Additionally, emission is still observed at the same 
wavelength in Fm(A). When Fm(A) (redward of Hy) is fitted 

6 Because narrow emission -possibly associated with Hy- is also 
present at this wavelength in F M (A) and F M (C), it might be that only 
a fraction of [O m] ,14363 is micro-lensed. 
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Fig. 9. Zoom on the 3080 < A < 3540 A range (rest frame) for 
^M/j(A) (dashed red) and Fmh(C) (thin solid blue) and Fm(A) 
(thick solid red) and Fm(C) (dashed-dotted blue). The pseudo 
continuum model of Fm^ is overplotted in short dashed pink 
lines. The emission in Fm(A) is not compatible with and is 
likely Fe n emission. However, there is a conspicuous excess of 
flux in F Mfl {A) w.r.t. the model (dotted pink line) which might 
also be due to Fe n. For image C, all the emission in that range 
is micro-lensed. 

with the Fe ii opt template, we find a lack of emission w.r.t. the 
template in the range 4630-4800 A. This suggests that the ex- 
cess of flux in Fm[i{A) corresponds to Fe n emission missing in 
Fm(A). In the range 4630-4800 A, Fen is mainly produced by 
the four following multiplets: [Fe n] 4F, [Fe n] 20F, Fe n 50 and 
Fen 43 (Veron-Cetty et al. 2004). If the micro-lensed Fen is 
associated with a particular Fe n transition, one would expect 
to observe other components of the multiplet associated with 
that transition elsewhere in F Mtl (A). Among the 4 mentioned 
multiplets, Fen 50 and Fen 43 only emit in the range 4630- 
4800 A, but [Fe u] 4F is also emitting at A 4889 A and [Fe n] 20F 
is also emitting at AA4874, 4905 A and around [O m] A/14959, 
5007. Since these other components of the multiplet are sig- 
nificantly blended with the [O m] and H ft emissions, we indeed 
cannot identify which are the micro-lensed Fe n transition(s). In 
conclusion, the present data indicate that at least a fraction of 
Fe n in the range 4630-4800 A is emitted in a compact region, 
possibly as compact as the VBC of Mg n. 

Similarly to what is observed in the range 4630-4800 A, an 
excess of flux is also observed in the rest-frame range 3080- 
3540 A (where no Fen templates exist). Indeed, we can see 
in Fig. [9] that a fraction of the emission in that range is mi- 
crolensed in A & C while another fraction is micro-lensed only 
in image C. Although it is known that the emission in that 
range is mainly due to Feu, we cannot demonstrate that the 
micro-lensing behaviour we observe is associated with two dif- 
ferent Fen emitting regions (as for the range 4630-4800 A) or 
not. If not associated with Fe n, we have no credible emission 
candidate 7 for the fraction of the emission (in the range 3080- 
3540 A) micro-lensed in A & C. This open question needs care- 

7 Due to the steep flux increase observed in the range ~ 3040- 
3 120 A, the 3080-3540 A range was imperfectly modelled with our 



ful identification of the Fe n emission -and absorption- lines in 
the range 3080-3540 A. This is beyond the scope of the present 
paper. 

7. Additional results 

7.1. The NEL extension 

The spectrum extraction of A, B and C has revealed that the 
spatial profile of the [Om] A/14959, 5007 region is significantly 
different from the one of the underlying continuum especially 
for images A and C. Indeed, after extraction, significant resid- 
uals are present around both images A & C while nothing com- 
parable is observed around image B (Sect. 1231 ). This indicates 
that the flux in the NEL is partly resolved along the slit di- 
rection. On the other hand, we did not find a significant offset 
(i.e. offset < 1 pixel) of the centroids of the NEL images w.r.t. 
to the centroid of the continuum emission. Using a Singular 
Isothermal Ellipsoid + external shear (SIE+y) to model the 
lensing galaxy (Paper I), we calculate that the absence of such 
an offset implies an upper limit of roughly ~ 30/T 1 pc on the 
off-centering of the [O m] emission w.r.t. the continuum emis- 
sion. Once we have noticed the absence of an off-centering 
of the [O in] /L14959, 5007 emission, we can safely impose a 
limit on its size. Indeed, using the SIE+y lens model, we find 
that images A, B and C are merging for a source size radius 8 
> 1 lO/i 1 pc. This size is compatible with the typical size of 
the NLR (~ 1 kpc; Bennert et al. 120021 1 . Stronger constraints 
could be derived by combining integral field spectroscopy of 
J 1 1 3 1 with a proper method to estimate uncertainties due to the 
seeing, sampling rate and lens model degeneracies (Yonehara 
120061 

7.2. The intrinsic absorption lines 

Another remarkable result is the presence in our spectra of a 
Mg n absorption doublet blueshifted at z =0.654 . According 

to its velocity shift (Av 660 km s~'), this absorption system 

can be classified as an "associated absorption line" (AALs; see 
e.g. Hamann & Sabra l2003l) . 

These lines disappear from the spectrum ratios A/B and 
C/B. This indicates that the absorbed flux is proportional to the 
flux coming from the continuum+BLR. This implies that the 
region at the origin of the absorption must cover both the con- 
tinuum and the BLRs and that, within the uncertainties, their 
depths are identical in the spectra of the three images A, B, and 
C. This is clearly illustrated in Fig. [10] where the absorption 
lines are normalized to a local effective "continuum" which in- 
cludes the broad Mg n emission line. 

The relative intensities of the blue and red lines of the dou- 
blet show evidence for partial coverage. Indeed if 4iue and 7 re d 
are the normalized residual intensities of the two lines of the 
doublet and r the optical depth, we expect /blue = e~ T and 

Balmer continuum model, suggesting that the microlensed flux in 
^mji(A) is not related to the Balmer continuum. 

8 We should notice that this estimate is model dependent. If e.g. 
multipoles are added to the lens model (Paper II), one finds that images 
A, B and C are merging for emission regions > 36/i~' pc. 
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Fig. 10. Associated Mg n absorption doublet normalized to a 
local effective continuum (i.e. continuum including the broad 
Mg ii emission) for the 3 lensed images A, B and C. The sim- 
ilar depth in the three images indicates that the origin of the 
absorption covers both the continuum and the BLR. 



Ired = e~° 5t in the case of full uniform coverage. The fact 
that < /blue as clearly seen in Fig. QJJ] suggests that only 
a part of the local continuum is covered by the absorption line 
region. Assuming homogeneous partial coverage (Hamann & 
Sabra 2003), we find from these intensities that the covering 
factor of the absorption region is around 20% (see e.g. the rela- 
tions given in Hutsemekers et al. 2004). This suggests that the 
absorption is intrinsic to the quasar, and exterior to the BLR. 

One can see the differential microlensing at work in Jl 131 
as a probe of the inhomogeneities in the absorbing medium. 
Indeed, by lensing more (less) strongly some regions of the 
source, differential micro-lensing increases (decreases) the 
contribution of a fraction of the intervening absorber to the total 
absorption. The nearly identical absorption depths seen in the 
three images (Fig. IToT > indicates that both the spatial distribu- 
tion and the optical depth of the absorbing clouds must be quite 
homogeneous over the continuum and BLRs. This is compati- 
ble with an absorption region composed of a large number of 
small absorbing clouds, their projected sizes being significantly 
smaller than the continuum region. A more quantitative study 
of this effect needs modelling of the absorbing medium, of the 
micro-lensing event as well as higher spectral resolution and 
S/N data. 



7.3. The lensing galaxy and neighbours 

The slit orientation #2 provides us with a medium resolu- 
tion spectrum of the lensing galaxy and of two neighbour 
galaxies gal #1 and #2 located at 55" and 95" from the lens 
(Sect. 12731 . We confirm the identification of the lensing galaxy 
by Sluse et al. (120031 ) as an elliptical galaxy at z = 0.295. 
Unfortunately, the S/N is too poor to derive the central veloc- 
ity dispersion. The galaxy gal#l is identified on deep /?-band 
imaging as a face-on spiral galaxy. There are 2 strong emis- 
sion features with an intensity ratio 2:1 observed at /M8461.5 
& 8488 A and a fainter one at A 8440 A. We identify these 
lines as Hor/t6562 and [Nn] /M6548, 6583 emission associated 
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Fig. 11. Smoothed spectra in the blue and red ranges (thin solid 
red) of the lensing galaxy (bottom panel), gal #1 (middle panel) 
and gal #2 (top panel) and associated redshifted template spec- 
tra (thick solid blue). A flux offset has been imposed for legi- 
bility. Due to mis-centering of gal #1 and #2 in the slit, there is 
significant slit losses for 5500 < A < 6500 A and A > 8500 A. 
Best redshifts are found to be z- 0.295 (Lens), z= 0.289 (gal 
#l)and z= 0.105 (gal #2). 



with intense star formation. A careful look at the smoothed 
spectrum of the galaxy reveals emission of [S n] AA6116, 6731, 
[On]/l3727 and Hy/14340 that definitely confirm a redshift 
value of z — 0.289 + 0.002 for this galaxy. Once smoothed with 
a Gaussian kernel, the spectrum of the second galaxy (gal#2) 
clearly reveals several absorption features among which Ca K 
and H AA 3933.7, 3968.5 A. The redshift of gal #2 derived from 
these absorption lines is z — 0.105+0.002. The spectra of these 
galaxies compared to the associated redshifted templates are 
displayed in Fig.fTTI 

Gal#l and gal#2 are in the direction of the external shear 
(with the convention that the shear y points towards or oppo- 
site to the perturber). If we model the perturbing galaxy as 
a Singular Isothermal Sphere, we can calculate the amount 
of shear caused by that perturber (e.g. formula A. 20. of 
Momcheva et al. 120061 ). Using a conservatively large value <x = 
300 km s _1 for the velocity dispersion, we estimate that the 
amount of shear caused by gal#l and gal#2 is only y ~ 0.01 1 
for each galaxy 9 . Even then, this is not sufficient to explain the 
shear value y ~0.1 necessary in the lens models (e.g. Papers I 
& II). This suggests that these two perturbers might be mem- 
bers of more massive galaxy groups. Interestingly, Williams et 
al. (12006b -based on two band imaging and red sequence find- 
ing technique- found evidence for two galaxy groups likely at 
z =0. 1 9 and z =0.29 in the field of J 1 1 3 1 . The redshift of gal#2 
together with the measured redshift z - 0.1035 and z = 0.1006 
of two other field galaxies obtained during the Las Campanas 
redshift survey (Shectman et al. 1996 ) suggest a redshift z ~ 0.1 
for the foreground galaxy group. These two possible groups 



9 Although the angular separation from the lens is larger for gal#2, 
its shear amplitude is as high as for gal#l because of its smaller red- 
shift 
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might play a role in the modelling of that system. This will be 
further investigated elsewhere. 

8. Caveats 

In this section, we first further discuss the use of [O m] as a rea- 
sonable estimate of the macro-lens model flux ratios. Second 
we estimate whether the results might be affected by the time 
delay between A-B-C. Third, we briefly address the question 
of the effect of intrinsic quasar absorber that might modify the 
BEL profile independently of any micro-lensing. 

8.1. Macro-amplification ratios 

Although it seems reasonable to use the flux ratios in [O m] as 
an estimate of the macro-model flux ratio, one may ask how 
robust our conclusions are with respect to that assumption. 
For this purpose, we can compare the flux ratios in [Om] to 
those expected by a simple smooth lens model. Using the 
SIE+y model derived in Paper I as a fiducial model, we find 
Fa/Fb = 1-65 and Fc/Fg - 0.9 while the values measured 
in [Om] are F A /F B = 1.97 and F C /F B = 1.33. Although there 
are significant differences likely due to an imperfect modelling 
of the lens potential, the values measured in [Om] are more 
similar to the model than those measured in the BEL or in 
the continuum. Therefore, using the fiducial model estimates 
of the macro model flux ratios will not change our conclusions. 
Indeed, using the macro model values instead of the [O m] flux 
ratios as a proxy of the macro-amplification ratios would only 
change the amplitude of micro-lensing in A and C but nothing 
else. 

One may ask whether the [Om] emitting region is large 
enough to be unaffected by micro-lensing. Indeed, using equa- 
tion [9] one find that a region larger than lOOpc can be signif- 
icantly microlensed by substructures more massive than typi- 
cally 5 1O 7 M . Such a massive substructure is unlikely to af- 
fect image A. Indeed, if present, it should significantly mod- 
ify the BEL while we have evidence that the latter is only 
marginally de-amplified compared to the continuum. On the 
other hand, milli-lensing of image C is more likely since we 
observe a significant de-amplification of a large fraction of the 
BLR. Consequently, we may hypothesize that the NLR in C is 
partly demagnified. The only consequence of such an assump- 
tion would be that the effect of micro-lensing of the continuum 
and the BELs in C are stronger than estimated. Alternatively, 
magnification (instead of de -magnification) of the NLR seems 
incompatible with the data. Indeed, if the NLR is magnified by 
substructures, one expects that the BLR and the continuum will 
be more strongly affected, which is not observed. Thus, only an 
ad-hoc scenario in which the NLR is magnified by a massive 
substructure (e.g. a dwarf galaxy) and the BELs and continuum 
are de-magnified by stars (e.g. the stars in the dwarf) would be 
compatible with the data. 

8.2. Bias due to the time delay 

Recently, Morgan et al.(2006) argued that the time-delay be- 
tween A and B (C) is around 15 days. One may thus imag- 



ine that time variability between A and B (C) may play a role 
in the interpretation of our measurements. However, over the 
nearly 2 years of monitoring performed by Morgan et al., in- 
trinsic variability never exceeded 0.05 mag over a period of 
15 days. Intrinsic variability may thus be considered as an ad- 
ditional source of noise affecting our measurements at a level 
of typically 5%. This is not only true for the continuum but 
also for the BEL that is mainly formed through photoionisation 
processes which imply that the BEL flux variations respond 
linearly to the flux variations of the continuum (e.g. Peterson et 
al. [1985). Consequently, the time-delay between A and B (C) 
should only marginally affect our results. 

8.3. Bias due to intrinsic absorbers 

Green (2006) recently investigated how the small line of sight 
differences existing between lensed images may modify the ob- 
served emission line profiles. Green suggests that some warm 
absorber outflow located close to the QSO continuum might 
modify the BEL line profile since it would be present along 
the line of sight of one image but not of the other. However, if 
intrinsic absorbers modify the BEL in A-B-C it is unrealistic 
that they would affect the whole emission profile. Partial "de- 
formation" of the broad line profiles due to a warm absorber 
might bias the Ffi decomposition of the BELs (Sect. I3.21 l. be- 
cause the method assumes a similar intrinsic profile in each 
component. The similarity between the B and C line profiles 
suggests that BEL differences due to line-of-sight effects are 
not significant in these images. If present, this effect likely oc- 
curs in image A for which the H/3 line profile appears more 
asymmetric. However, H/3 absorbers are very rare: up to now, 
broad absorption lines in H/3 have been observed only in 4 
systems (Hall 2006). Long term spectroscopic monitoring will 
shed light on such an effect if present at all. 

9. Conclusions 

This paper is devoted to a thorough analysis of the long slit 
spectra obtained for the gravitationally lensed quasar Jl 131- 
1231. The spectra of the three brightest images A-B-C en- 
abled us to estimate the flux ratios in different emitting regions 
(namely the continuum, the broad line and the narrow line 
emitting regions) and unveil the micro-lensing effects occur- 
ring in that system. Evidence for differential micro-lensing of 
the broad emission lines enabled us to perform a phenomeno- 
logical study of the structure of the Broad Line emitting Region 
(BLR). 

Due to the better resolving power of the present spec- 
tra, new redshift estimates of the source and of the lensing 
galaxy have been performed. Based on the Mg n emission, a 
systemic redshift z =0.657+0.001 has been deduced. The nar- 
row emission lines (except [O u] /13727; z = 0.656 + 0.001) are 
blueshifted w.r.t. the systemic redshift. A redshift z = 0.654 + 
0.001 is measured for those lines. On the other hand, we con- 
firm the redshift of the lensing galaxy to be z = 0.295 + 0.001. 
Finally, we have estimated the redshifts of the 2 galaxies gal #1 
and #2 located at resp. 55" and 95" from the lens to be resp. 
0.289 + 0.002 and 0.105 + 0.002. 
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The results derived from the micro-lensing analysis and its 
implications for the quasar structure are summarized below: 

1. Different flux ratios were derived for the continuum, the 
Broad Emission Lines (BELs) and the Narrow Emission 
Lines (NELs), indicating that micro-lensing is at work in 
at least 2 images. The simplest scenario explaining our ob- 
servations consists in micro-lensing de-amplification of im- 
ages A and C. Additionally, we found evidence for differ- 
ential micro-lensing of the BLR I0 . Since a larger fraction 
of the BELs is micro-lensed in image C, the Einstein radius 
of the micro-lens is likely larger in C than in A. Using the 
relation between the size of the BLR and the luminosity of 
the QSO (Kaspi et al. [20051 Bentz et al. [2006b . we found 
that a micro-lens of a few solar masses is sufficient to affect 
the whole BLR. 

2. We found that the Mgn and the Balmer emission lines 
cannot be represented with the same velocity components. 
Especially, the very broad component of the Mg n emission 
is not present in the Balmer lines. Since it is micro-lensed 
in image C and A, this component likely comes from a very 
compact region. 

3. The differential micro-lensing of the BELs confirms that 
the size of an emission line region is anti-correlated with 
the FWHM of the corresponding line component. 

4. We have an indication that differential micro-lensing in im- 
age A resolves the BLR in velocity. Monitoring of such an 
effect combined with modelling of the BLR offers interest- 
ing avenues to probe the BLR structure and kinematic. 

5. We found that the narrow emissions observed in the Balmer 
lines have the same characteristics as the [O m] AA4959, 
5007 narrow emission, namely same width, redshift and 
magnification ratio. This argues in favour of a common 
emitting region. On the other hand, we rule out the exis- 
tence of a similar narrow emission for Mg n. 

6. An intriguing narrow micro-lensed emission incidentally 
coincides with the [O m] /14363 emission, suggesting a very 
compact emission region for that line. The robustness of 
this result is however questionable. Indeed, this emission 
might also be produced by the differential micro-lensing 
of H j. Additionally, one would expect to see a similar be- 
haviour for e.g. [Ne v] emission, which is not observed. This 
result clearly needs further investigation. 

7. The [O m] /1/14959, 5007 narrow emission lines are partly 
spatially resolved. Neglecting seeing, resolution effects, 
and lens model degeneracies, we argue that it provides us 
with a lower limit on the size of the NLR of 1 IO/1 1 pc, fully 
consistent with typical NLR size. 

8. We have firm evidence that a large fraction of the near-UV 
and optical Fe n emission arises in the outer parts of the 
BLR. We have however shown that a smaller fraction of the 

10 During the referee process, Sugai et al. (2007 ) published an anal- 
ysis of IFS data of Jl 131 in the HyS-[Om] -M4959, 5007 range. 
Although obtained at a different epoch (Feb. 2005), these data dis- 
play emission line flux ratios quite similar to those presented here and 
interpreted as micro-lensing of images A and C. Similarly to our re- 
sults, these data also suggest differential micro-lensing of the BLR as 
well as spatially resolved [O m] emission. 



Fe ii is emitted in the inner parts of BLR. Compact emission 
from Fe n is very likely identified in the rest frame range 
4630-4800 A. Micro-lensed pseudo continuum emission is 
present in the range 3080-3540 A. This emission is likely 
associated with Fe n. 
9. An associated absorption line doublet was observed in 
Mg ii at the same redshift as the NELs. We have found that 
the absorbing medium is intrinsic to the QSO and covers 
both the continuum and BLR emission, with a covering fac- 
tor of ~20%. Additionally, using differential micro-lensing 
as a probe of the inhomogeneities in the absorbing medium, 
we have shown that both the spatial distribution and the op- 
tical depth of the absorbing clouds are homogeneous over 
the BLR and the continuum. 

Although our phenomenological description of the mi- 
crolensing of the BELs is imperfect in several aspects, it nicely 
illustrates the wealth of information one can retrieve from spec- 
troscopic observations of micro-lensed QSOs. Clearly, more 
detailed models as well as coordinated spectroscopic (integral 
field) monitoring of lensed quasars where micro-lensing takes 
place will improve our understanding of both the micro-lensing 
phenomenon and of the QSO structure. 
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